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A search for supersymmetry is presented using a sample of events with b jets and 
missing transverse momentum. The search uses a data sample of proton-proton col- 
lisions at a centre-of-mass energy of 7 TeV, corresponding to an integrated luminosity 
of 35 pb"\ collected with the CMS detector. A total of 0.33 +033 (stat.) ± 0.13 (syst.) 
events is predicted, using control samples in the data, to arise from standard model 
processes, and one event is observed in the data. Upper limits are set at the 95% 
confidence level on the cross sections of benchmark supersymmetric models. 
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1 Introduction 

Supersymmetry (SUSY) |IlH5| is an extension of the standard model (SM) of particle physics, 
which can solve the "hierarchy problem" miZl and provide a candidate for cold dark matter |l8|. 
For a large class of supersymmetric parameter sets, squarks (q), the SUSY partners of quarks, 
are relatively light. In this case, significant event yields at the Large Hadron Collider (LHC) 
can result from strong production of squarks, which subsequently decay giving a weakly inter- 
acting lightest supersymmetric particle (LSP). If bottom and top squarks, which can decay to b 
quarks, are relatively light, there may be an abundance of events with one or more b-quark jets 
and momentum imbalance transverse to the beam line due to the undetectable LSPs. 

This Letter describes a search for events with two or more hadronic jets, at least one of which 
must be b tagged [9J, and significant transverse momentum imbalance. It extends a similar 
search without a b-tag requirement llTOl . The momentum imbalance is characterized |ITT| by the 
ratio of the pj of the second-highest-px jet and the invariant mass formed from the two highest- 
px jets. This ratio can be estimated by uj = \ ^J^^^\ , where = | E; Pj^' 1/ = Ei Pj^' > 

and pj ' is the momentum transverse to the beam line for jet / in an event. The jets in an event 

are grouped into two pseudo-jets and AHj is the minimal value of | — p^^'^"'^"^'^'^ | over 

all combinations; this approach optimizes rejection of backgrounds with apparent Ifj from 
instrumental effects and other sources. 

The main backgrounds are due to standard model multijet production (hereafter denoted "QCD 
background"), electroweak W and Z boson production (EWK), and top quark pair production 
(tt). Owing to low average Uj, the QCD background is effectively rejected by a requirement on 
Oij. The b-tag requirement further suppresses the QCD and EWK backgrounds. 

The results of the search are characterized in terms of the mSUGRA/CMSSM |IT2l[T3l scenario 
of SUSY. These models are described by four parameters and one sign: the universal scalar 
and gaugino mass parameters, mo and mi/2, respectively; the universal trilinear coupling, Aq; 
the ratio of the two Higgs doublet vacuum expectation values, tan /3; and the sign of the Higgs 
mixing parameter, sign(/^). Three signal points are considered as benchmarks: LMO, LMl, both 
discussed in Ref. [10], and LMB (corresponding to ntQ = 400 GeV, mi/2 = 200 GeV, Aq = 
GeV, tan/5 = 50, and sign(^) > 0), chosen to be near the edge of sensitivity of this search in 
mSUGRA/CMSSM parameter space. 

The analysis presented here uses a data sample of proton-proton collisions at 7 TeV, correspond- 
ing to an integrated luminosity of 35 pb^^, collected with the Compact Muon Solenoid (CMS) 
detector, at the LHC. The main components of CMS are a silicon pixel and strip tracker, the 
crystal electromagnetic calorimeter, and the brass/ scintillator hadron calorimeter, all placed in 
a 3.8 T axial magnetic field, complemented by gas-ionization detectors embedded in the steel 
return yoke, to measure muons. A detailed description of the detector and its performance can 
be found in Ref. [14] . In the cylindrical coordinate system of CMS, (p is the azimuthal angle 
and the pseudo-rapidity {?]) is defined as // = — In [tan (0/2)], where 6 is the polar angle with 
respect to the counterclockwise beam direction. 



2 Event Selection 

The event selection requirements are mostly identical to those in Ref. [10|. Events in the search 
sample are collected with triggers based on Hj computed from jets reconstructed at trigger 
level. A muon-enriched control sample is collected with triggers requiring a muon. Events 
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3 Background Estimation 



must have a good reconstructed pp collision vertex [15]. Jets are reconstructed as clusters of 
energy in the calorimeters by the anti-kj algorithm [16J with a distance parameter of 0.5, and 
are required to have energy transverse to the beam, Ej, in excess of 50 GeV and \ri\ less than 3. 

To perform a fully hadronic final state search and to reduce the backgrounds, events with an 
isolated lepton (electron or muon) or photon are vetoed, and events consistent with having 
apparent Wj llTOl are rejected. Selected events are required to have at least two jets, both with 
Ej > 100 GeV, \rj \ < 2.5 for the highest-Ej jet, Hj > 350 GeV, at least one jet tagged as originating 
from a b quark, and ctj > 0.55. 

Jets are b tagged using a discriminator based on the impact parameter significance of tracks in 
a jet (Track Counting High Purity discriminator, TCHP |9|), with a "tight" selection (TCHP> 
3.41) designed to have a light-flavour contamination of less than 0.1%. Looser b-tagging se- 
lections are used to produce various control samples. An event is said to be anti-tagged if it 
contains no jets with a loose b tag (TCHP > 1.19). To remain within the acceptance of the pixel 
tracker, only jets with a central axis of | ?^ | < 2.4 are considered for b tagging. 



3 Background Estimation 

The backgrounds for this search can be categorized into three main groups: namely QCD, 
EWK, and tt. The contamination from tt is mainly in the tau decay mode. The vast majority of 
events from the QCD background do not feature large transverse momentum imbalance and 
are therefore rejected by the (Xj > 0.55 requirement. The EWK backgrounds consist of W and 
Z boson production, with genuine missing energy due to decay neutrinos. The requirement of 
at least one b jet greatly reduces the EWK and QCD backgrounds. The dominant background 
for the analysis arises from tt production, in which b jets and genuine missing energy due to 
neutrinos can arise from the top quark decay chains. 



A procedure based on control data samples, described in Section 3.1 is employed to estimate all 
backgrounds simultaneously. In this method, the fraction of all events with ocj > 0.55, denoted 
F{aj > 0.55), is measured in a lower-Hx control region and applied in the signal region. 



The Z — > vt7 and tt background yields are cross-checked separately, as discussed in Section 3.2 
The tt cross-check uses muons to emulate the hadronic decays of taus. The cross-check of 
Z — > 1/1/ utilizes Z — > events for which dli is determined after excluding the muons. 

3.1 Background Prediction Using cci vs Hj Extrapolation 

In SM simulation studies [10], V{cii > 0.55) has no Hj dependence in events with large genuine 
missing transverse energy, i.e., the tt and EWK backgrounds. In the QCD background, how- 
ever, F(ax > 0.55) is expected to be a decreasing function of Hj because of the Hj dependence 
of the factors contributing to apparent Uj, such as jet energy resolution and jet Ej threshold 
effects. 

In data control samples, F(ft:x > 0.55) is consistent with having no Hj dependence, which 
indicates that the tt and EWK backgrounds dominate. The larger anti-tagged data sample 
is also consistent with having no Hj dependence. Because a tight b-tag requirement further 
suppresses the QCD background, the tight tagged data sample is expected to have a negligible 
QCD contribution and therefore ¥{ixj > 0.55) independent of Hj. 

The total background is estimated by measuring F(aT > 0.55) = 1.48^^ 43 ^ ^ control 

region with 250 < Hj < 350 GeV and multiplying this fraction by the number of events in the 
signal region before the aj > 0.55 requirement. In data, this procedure yields a prediction of 



3.2 Cross-Checks of Z ^> yy and tt Background Contributions 
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Table 1: Predicted and observed numbers of events for 35 pb ^. The prediction comes from 
the Dij vs Hj extrapolation described in Section [ 



3.1 



N-jets 
> 2 



Background Prediction Data 
0.33 toil (stat.) ± 0.13 (syst.) 1 



LMO 
14 



LMl 
2 



LMB 

5 
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Figure 1: The ocj distributions for Z 
emulation of hadronic tau decays (dashed red) 



emulation of Z — > (solid blue) and muon 



0.33 lg33 (stat.) ±0.13 (syst.) events. The statistical uncertainty is dominated by the presence of 
one event with ocj > 0.55 in the control sample. The systematic uncertainty on the prediction is 
given by the difference in F(aT > 0.55) measured in the tight and loose tagged control samples. 
Table [ijlists this background prediction, the observation in data, and the expected contribution 
of SUSY signal for points LMO, LMl, and LMB. 

3.2 Cross-Checks ofZ^vv and tt Background Contributions 

While the above background estimate is the one used in this search, we perform auxiliary 
measurements to cross-check the Z — > and tt background components, which together are 
expected to comprise the majority of the background. As would be crucial in case of an ob- 
served excess, these cross-checks provide an overestimate of the Z — > vv and tt background 
components. 

For Z — ?► vv, a sample of Z — ?► events is selected with two or more jets but no dlj, Hj, or b- 

tagging requirements. The solid blue line in Figure [ijshows the oci distribution for the resulting 
events. The fraction of these events containing a b-tagged jet is measured. Then, a sample is 
selected with no b-tag requirement, jet Ej = 75 GeV thresholds on the two highest- jets, Hj 
> 275 GeV, and oci > 0.52. The number of events in this sample is scaled by the measured b-tag 
fraction in the other sample, corrected for the muon identification efficiency and acceptance, 

and multiplied by the ratio of branching fractions — > vv)^ ^ ^ This procedure gives 

BR(Z — y ji n ) 

an overestimate of the number of Z — )■ i/v events in the signal region owing to less stringent 
requirements than in the final selection, and yields 0.48 ± 0.39 events. 

Simulation studies indicate that most of the tt background comes from events with hadronic 
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4 Signal Selection Efficiency 



Table 2: Cumulative and individual efficiencies for the selection in three SUSY benchmark 
points. For each point, the left and right columns represent the individual and cumulative 
efficiencies, respectively. Different benchmarks have different b-tag efficiencies due to different 
average numbers of b quarks per event. The fraction of events containing at least one b quark 
before the b-tag selection is 66% in LMO, 18% in LMl, and 91% in LMB. 



Requirement 


LMO 


LMl 


LMB 


Pre-selection 


98% 


98% 


98% 


98% 


98% 


98% 


Lepton/Photon Veto 


57% 


56% 


55% 


54% 


61% 


60% 


Jet Requirements 


51% 


28% 


63% 


34% 


54% 


33% 


Ht > 350 GeV 


90% 


25% 


94% 


32% 


97% 


32% 


Trigger 


99% 


25% 


99% 


32% 


99% 


31% 


Apparent Uj Veto 


68% 


17% 


81% 


26% 


65% 


20% 


Tight b-tag 


31% 


5.3% 


12% 


3.0% 


54% 


11% 


aj > 0.55 


14% 


0.7% 


29% 


0.9% 


12% 


1.3% 



tau decays. To estimate the hadronic tau decay yield, F{ocj > 0.55) is first measured in a sample 
with Ej = 80 GeV thresholds on the two leading jets, Hj > 280 GeV, at least one medium b- 
tagged jet (TCHP > 1.91), and one or two muons. These selection requirements are chosen to 
be less strict than the signal selection in order to increase the number of events in this sample. 
The muons are used to emulate the hadronic decays of taus. To do so, for each muon the 
presence of a tau jet is emulated with an Ej value set to a fraction of the muon pj, using a 
distribution taken from simulation. The dashed red line in Figure [l] displays the resulting ocj 
distribution. The measured value of F(aT > 0.55) in this sample is multiplied by the number 
of emulated events in the signal region before the aj requirement. This value is corrected for 
the muon selection efficiency, acceptance and the hadronic tau decay branching ratio to obtain 
the hadronic tau decay yield. The predicted hadronic tau decay yield is increased by 38%, 
as determined in simulation, in order to account for the entire tt background. The procedure 
yields a 25% overestimate of the total tt background in simulation. In data, 1.4 ± 0.5 events are 
predicted. 

4 Signal Selection Efficiency 

To interpret the results of this search in terms of a given signal model, the selection efficiency for 
that model must be determined. Table |2]lists the cumulative and individual efficiencies for the 
event selection in the three SUSY benchmark models LMO, LMl, and LMB, from which events 
are generated at leading order (LO) via pythia 6.4, tune Z2 |1T7| using parton distribution 
functions provided by CTEQ6.6 [18J. Without b tagging, the cumulative efficiencies for LMO 
and LMl are about 85% of those in Ref. [lOJ, because of a more stringent lepton and photon 
veto. Table|3]lists the relative systematic uncertainties on the signal yield, which are dominated 
by the uncertainty on the b-tagging efficiency, described below. The other uncertainties and the 
methods used to obtain them are similar to Ref. IITOl . 

The b-tagging efficiency is measured from inclusive dijet events in which one jet has an asso- 
ciated muon and another "away" jet has a TCHP value of at least 1.0. The relative fraction of 
jets from b quarks in a data sample is determined by a fit to the distribution of transverse mo- 
mentum of muons relative to their associated jet axis, Elllll/ which is larger for jets from 
b quarks than from other flavours. This fit is to a linear combination of simulation-derived 
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Table 3: Relative systematic uncertainties on the signal yield. 



Source 


Uncertainty (%) 


Luminosity 


4 


JES 


3.5 


Jet Energy Resolution 


1 


Trigger Efficiency 


1 


Apparent Ut Veto 


4 


Lepton/Photon Veto 


4 


b-tag Efficiency (LMB) 


20 


Total 


22 



Pj expected distributions from different flavours. The fitted b fractions for jets passing and 
failing the analysis b-tagging requirement are used in the b-tagging efficiency calculation. This 
efficiency is measured separately for jets with \ f]\ > 1.4 and \ f]\ < 1.4, in four ranges of jet £x- 
The ratio between the b-tagging efficiency measured in data and in simulation is taken as the 
efficiency scale factor for a particular range in Ej and \ f]\. 

Systematic uncertainties on the scale factors arise from potential biases in the ' fitting pro- 
cedure. These uncertainties are measured by varying the muon-to-jet matching and muon pj 
thresholds, fraction of gluon splitting to bb, jet energy scale and resolution, jet angular reso- 
lution, and b-tagging requirement on the away jet. The effect of measuring the scale factors 
using only semi-leptonic b decays is also accounted for. The scale factors are used to correct the 
expected event yield at each signal point for differences between the efficiencies in data and 
simulation. For example, for LMB the application of the scale factors translates into a change in 
the yield by a factor 0.87 ± 0.18. The systematic and statistical uncertainties give a total relative 
uncertainty of 20% in LMB, with a similar uncertainty of 23% in LMl. 



5 Results 

The observation of one data event in the signal region is consistent with background expecta- 



tions. Combining the expected signal and background prediction from Section 3.1 and using 
frequentist statistical methods in the manner of Ref. II20I with the Profile Likelihood ratio l|2T| 
to handle nuisance parameters, we derive 95% confidence level (CL) cross-section upper limits 
(tT-*^) of 18.9, 15.4, and 10.2 pb for LMO, LMl, and LMB, respectively The effect of possibly 
overestimating the background due to signal contamination in the control regions increases 
the (7^1^ value to 22.1 pb for LMO, 16.7 pb for LMl, but is negligible for LMB. To quantify the 
sensitivity with reduced dependence on the amount of b-quark production, a 95% CL upper 
limit on the cross section times branching ratio to at least one b quark of 4.0 pb is determined 
in LMl. 

The resulting excluded region in the (mo, mi/2) plane for a reference model with CMSSM pa- 
rameters Aq = GeV, tan /3 = 50, and f/ > is shown in Figure |2] The expected and ob- 
served exclusion regions are calculated using next-to-leading-order (NLO) cross sections, ob- 
tained with the program Prospino [22 J. The excluded region is extended with respect to that 
of Ref. [lOJ without b tagging, also shown, for scenarios with increased b production, such as 
those with mo above 350 GeV. 
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6 Summary 



CMS Li„, = 35 pbWs = 7 TeV 
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Figure 2: Exclusion regions in the (mo, mi/2) plane for one set of CMSSM parameters, for this 
analysis (solid red), and the non-b tagged version fflOl (dashed green). 

6 Summary 

A search for events with multiple jets, at least one of which is b tagged, and significant trans- 
verse momentum imbalance has been presented. One event is observed, which is consistent 
with background expectations. The dominant background comes from tf production. The re- 
sults of the search are characterized as an exclusion region in CMSSM parameter space and 
95% CL upper limits on representative scenarios with expected cross section of approximately 
15 pb. The sensitivity of this search surpasses that of the Tevatron experiments 11231 124| and is 
comparable to a recent fully hadronic, b-tagged search from ATLAS I25j . 
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